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Abstract 



The dependence of the photon structure on the photon virtuality, Q 2 , is studied 
by measuring the reaction e + p — > e + + jet+ jet+ X at photon-proton centre-of-mass 
energies 134 < W < 223 GeV. Events have been selected in the Q 2 ranges « GeV 2 , 
0.1-0.55 GeV 2 , and 1.5-4.5 GeV 2 , having two jets with transverse energy E^ e * > 5.5 
GeV in the final state. The dijet cross section has been measured as a function of the 
fractional momentum of the photon participating in the hard process, x® . The 
ratio of the dijet cross section with x® BS < 0.75 to that with x® BS > 0.75 decreases 
as Q 2 increases. The data are compared with the predictions of NLO pQCD and 
leading-order Monte Carlo programs using various parton distribution functions of 
the photon. The measurements can be interpreted in terms of a resolved photon 
component that falls with Q 2 but remains present at values of Q 2 up to 4.5 GeV 2 . 
However, none of the models considered gives a good description of the data. 
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1 Introduction 



The photon at high virtuality, Q 2 , is commonly considered to be a point-like probe of the 
structure of a particular hadronic target []T|. However, the real photon (Q 2 ps GeV 2 ) 
has itself a partonic structure, which has been studied in two-photon reactions from e + e~ 
scattering @, and in jet production at HERA f|, ||. In this paper, the transition between 
the real photon and the virtual photon is investigated for < Q 2 < 4.5 GeV 2 using dijet 
events in ep scattering at HERA. 

The photon, in general, may have both a partonic structure and a point-like coupling 
to charged quarks and leptons. As a result, two types of process can contribute to jet 
production in ■yp interactions in leading order (LO) perturbative QCD (pQCD): the direct 
process, in which the photon couples directly to quarks at high transverse momenta, one 
of which scatters from a parton in the proton, and the resolved process, where a parton 
from the photon scatters from a parton in the proton. Conventionally, two types of 
resolved photon process are defined. In the first, the photon acts via an intermediate 
meson-like hadronic state whose description is essentially non-perturbative, so that a 
phenomenological parton density function must be introduced. In the second, the photon 
interacts initially by splitting into a qq pair at moderate transverse energy, a point-like 
perturbative process which is termed 'anomalous' and can in principle, for Q 2 > GeV 2 , 
be summed to all orders. The boundary between the two types of resolved process is 
factorisation-scale dependent. 

At a given photon virtuality, Q 2 , and hard QCD scale, /x 2 , both types of resolved process 
can in principle occur. It is usually accepted that, at low Q 2 , the hadronic type is impor- 
tant, while at higher Q 2 , resolved processes are dominated by the anomalous type. The 
general expectation is that the contribution to the dijet cross section from both types of 
resolved photon processes should decrease relative to the contribution from direct photon 
processes as the virtuality of the photon increases towards /i 2 , i.e. the partonic content 
of the photon becomes suppressed [§, f| 0, §■ The first measurements came from the 
PLUTO collaboration M. The HI collaboration has also studied the transition between 
photoproduction and deep inelastic scattering by measuring, in the 7*p CM frame, in- 
clusive jet cross sections for real and virtual photons JTU] and dijet cross sections |TT| for 
Q 2 > 1.6 GeV 2 . 

The resolved and direct components can be separated on the basis of the variable x® BS , 
which is the fractional momentum of the photon partaking in the production of the dijet 
system. This variable is defined as: 

OBS _ 2^ jets e 



2yE e 

where and rf et are the transverse energy and pseudorapidity of the jet defined in 
the laboratory frameQ. The variable y is defined as y = 1 — 2^r(l — cos# e ), where E e is 

1 The ZEUS right-handed coordinate system is defined with the origin at the nominal interaction point 
by the Z axis pointing in the proton beam direction and the X axis pointing horizontally towards the 
centre of HERA. 
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the positron beam energy and E' e , 0' e are the energy and polar angle, respectively, of the 
scattered positron. Since x° BS is well defined at all orders in pQCD, measurements based 
on it can be compared with theoretical predictions at any given order. 

At x^ BS > 0.75, the direct component dominates, while the x® BS < 0.75 region is sensitive 
mainly to the resolved component. However, events with low values of x® BS can also be 
produced when initial- and final-state parton showers give rise to hadronic activity outside 
the dijet system. 

Experimental x° BS distributions obtained from 1995 ZEUS data are presented in this 
paper. The ratio of the measured cross sections for x^ BS < 0.75 and > 0.75 is then given 
as a function of Q 2 . The values of the ratio are compared with theoretical calculations 
at both LO and next-to-leading-order (NLO) pQCD computed using the JetViP program 



2 Experimental Setup and Data Selection 

During 1995, HERA operated with protons of energy E p = 820 GeV and positrons of 



energy E e = 27.5 GeV. The ZEUS detector is described in detail elsewhere |L3], |i4f . 
The main components used in the present analysis are the uranium-scintillator sampling 
calorimeter (CAL) f!5 |, the beam pipe calorimeter (BPC) ||16|| , and the central tracking 



chamber [T7| positioned in a 1.43 T solenoidal magnetic field. The CAL energy resolution 



for positrons, under test beam conditions, was measured to be 0.18/y -Eg(GeV). The point 
of impact of the positron in CAL was measured with a resolution of 3 mm, resulting in a 
Q 2 resolution of 8%. The systematic uncertainty on the absolute value of E e is 2%. The 
BPC was installed 294 cm from the interaction point in the positron direction in order 
to tag scattered positrons at small angles (15-34 mrad). It measured both the energy, 
E e , of the scattered positron and the position of its impact point. The energy resolution 

of the BPC is 0.17/ J E' e (GeV) and the position resolution is 0.5 mm, resulting in a Q 2 



resolution of 6%. The systematic uncertainty on the absolute value of E' e is 0.5%. 



The events were selected online via a three-level trigger system |14], [19] using the same 
selection algorithms as in a previous dijet publication ||], except that in the third-level 
trigger (TLT) the events were required to have at least two jets with E^ LT > 4.0 GeV 
and rj TLT < 2.5. The sample was separated offline into subsamples corresponding to three 
different Q 2 ranges: 

• Events with quasi-real photons (Q 2 ~ GeV 2 , named PHP in the following) were 
selected by requiring that no identified positron was found in the CAL with energy 



E e > 5 GeV |2(J and y < 0.7. The resulting sample had Q 2 < 1.0 GeV 2 with an 
estimated median of 10~ 3 GeV 2 ; 

Events at intermediate Q 2 (IQS) were selected by requiring that the scattered 
positron was measured by the BPC. In this data set, the BPC tagged events with 
photon virtualities in the range 0.1 < Q 2 < 0.55 GeV 2 . For this sample the energy 
of the scattered positron was required to be E e > 12.5 GeV; 



2 



Deep inelastic scattering (DIS) events at low Q 2 (LDIS) were selected by requiring 
that the outgoing positron was measured in the CAL. The energy of the scattered 
positron was required to satisfy E e > 11.0 GeV |20| and the Q 2 range was restricted 
to 1.5 < Q 2 < 4.5 GeV 2 . 



For all three samples, additional cuts of 0.15 < yjs < 0.45 were applied, where yjs is an 
estimator^ of y [[Hj]. Due to the energy lost in the inactive material in front of the CAL 
and to particles lost in the rear beampipe, yjB systematically underestimates the true 
y by approximately 20%, an effect which is adequately reproduced in the Monte Carlo 
simulation of the detector. The combination of yjB and E e cuts ensured that all three 
samples corresponded to the same true y range (0.2 < y < 0.55). 

The longitudinally invariant kx algorithm |22| , in the mode described in a previous pub- 
lication [H, was then applied to the CAL cells to search for events with two jets in the 
final state. In the LDIS sample, the cells associated with the positron were excluded from 
the jet search. The two jets with the highest transverse energy were required to have 
pseudorapidity between —1.125 < rf^ < 2.2 and transverse energy Ef > 5.5 GeV. 

After all cuts, the jet search resulted in a sample of 58224 dijet events for the PHP 
sample, 353 dijet events for the IQS sample and 1172 dijet events for the LDIS sample. 
Approximately 10% of the events in each of the three samples had three or more jets. 
The PHP, IQS and LDIS samples correspond to integrated luminosities of 3.1, 3.3 and 
4.9 pb -1 , respectively. 



3 Data Corrections and Systematics 

The data were corrected for acceptance, smearing and kinematic cuts using the HER- 
WIG 5.9 |23| Monte Carlo (MC) model. Leading-order resolved (LO-RES) and direct 
(LO-DIR) processes were generated separately. Resolved photon events were generated 
not only in the PHP, but also in the IQS and LDIS regimes. The minimum transverse 
momentum of the partonic hard scatter (p™ m ) was set to 2.5 GeV. The GRV LO [^4] and 



the MRSA |35|| sets were used for the photon and proton parton distribution functions 
(PDF), respectively. To simulate possible interactions between the proton and photon 
remnants ('underlying event'), the option of multiparton interactions (MI) |2"E| , ^7] was 
included for the PHP sample. It has been shown that the simulation of the underlying 
event with MI improves the description of the energy flow around the jet axis for jet 
production from quasi-real photon-proton interactions [|J. 

The Monte Carlo events were processed through the full ZEUS detector simulation using 
the same cuts as applied to the data. The normalisations of the LO-RES and LO-DIR 
processes were extracted from the data using a two-parameter fit to the uncorrected x® BS 
distributions. This procedure was applied separately for each Q 2 range. 



2 Ujb — J2i(Ei ~ Ezi)/2E e , where Ezi — E{ cos6>; and Ei is the energy deposited in the CAL cell i 
which has a polar angle 9i with respect to the measured Z-vertex of the event. The sum runs over all 
CAL cells excluding those associated with a detected scattered positron. 
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Figure [l] shows uncorrected distributions of x° BS for PHP, IQS, and LDIS dijet events 
compared with the HERWIG simulation. Events both at high x° BS , associated mainly 
with direct photon processes, and at low x° BS , associated mainly with the resolved photon 
processes, are present in all Q 2 ranges. 

In the low x® region, the PHP data disagree with the simulation. Disagreement is also 
observed in the rf et , yjs, and E^ 1 distributions (not shown) and can be attributed to 
the presence of underlying event effects or uncertainty in the PDFs of the photon. The 
underlying event effects are most evident at low x® BS and low Ej?. For the combination 
of ]3™ m and photon PDF used here, the simulation of multiparton interactions does not 
reproduce the shape of the data in the low x° BS region ||. To take account of this 
disagreement in the correction of the data for migrations and acceptance, the Monte 
Carlo events have been reweighted as a function of x® at the hadron level so that the 
distribution agrees with the data. The result of the reweighting is shown in Fig. [1] (a). 
After the reweighting, the MC predictions for the ?7 jet , jjjb, and E^ 1 distributions also 
agree well with the data (not shown). 

The dijet differential cross sections, da/dx° BS , corrected to the hadron level, have been 
measured using the kx jet algorithm, in the three Q 2 regions with 0.2 < y < 0.55. The 
measurements have been made for two sets of jet transverse energy and pseudorapidity 
cuts: 



1. Low E J T et : E J T et > 5.5 GeV, -1.125 < ff et < 2.2 for both jets; 

2. High E 3 T et : E 3 T etl > 7.5 GeV, E 3 T et2 > 6.5 GeV, -1.125 < 7] jet < 1.875. 



The data with the low set of E^ 1 cuts are sensitive to the resolved photon component, 
but also to the effects of the underlying event. The data with the high Ej 1 cuts are not 
significantly influenced by the underlying event effects; this was established by means of a 
comparison (not shown) of the data with HERWIG without MI. The high E^ 1 cuts were 
chosen to be asymmetrical to facilitate a comparison with the NLO pQCD calculation. 

The cross sections at hadron level were obtained by applying a bin-by-bin correction to 
the measured dijet distributions binned in four x® BS bins (0.0625-0.25, 0.25-0.50, 0.50- 
0.75, 0.75-1.00) and four Q 2 bins (0.-1.0, 0.1-0.55, 1.5-3.0, 3.0-4.5 GeV 2 ). The correction 
factors take into account the efficiency of the trigger, the selection criteria and the purity 
and efficiency of the jet reconstruction. The efficiency and purity are determined as a 
function of i^ BS and Q 2 from the MC simulation [^] . In the PHP region, the correction 
factors lie between 1.25-1.43 and the purities between 0.50-0.64. In the IQS region, the 
correction factors are dominated by the BPC geometric acceptance and lie between 17.6- 
23.0 and the purities between 0.40-0.70. For the LDIS region the correction factors lie 
between 3.1-3.8, and the purities between 0.45 and 0.80. 

A detailed study of the systematic uncertainties of the measurements has been performed 



2^, |3(J. The uncertainties have been separated into those that are uncorrelated and 
therefore were added in quadrature to the statistical error and those that are correlated 
and presented separately. The uncorrelated systematic uncertainties originate from the 
residual uncertainties in the event simulation. The uncertainty associated with the E^ 1 
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cut is the dominant uncorrelated uncertainty for the PHP and IQS samples. When this cut 
is varied by the resolution of 14%, a systematic uncertainty between —9% and +12% 
results, except for x® BS < 0.25 in the IQS region where the uncertainty ranges between 
—28% and +10%. In the LDIS region, the dominant systematic uncertainty comes from 
the uncertainty in i^ BS (the x^ BS resolution is 0.05) and results in a systematic error 
between —2% and +6%, except for x° BS < 0.25 where the systematic uncertainty ranges 
between -25% and +36%. 

Two sources of correlated systematic uncertainties have been studied, one originating from 
the uncertainty of the CAL energy scale and the other from the use of different models for 
the description of the jet fragmentation process in the MC. The absolute energy scale of 



the jets in simulated events has been varied by ±5% |3l| . The effect of this variation on the 
dijet cross sections is ~ ±20%. The uncertainty associated with the jet fragmentation was 
studied by correcting the data to the hadron level using PYTHIA [[32| and comparing to 
the results obtained using HERWIG. The effect was estimated to be on average ~ 20%. 
In addition, there is an overall normalisation uncertainty of 1.5% from the luminosity 
determination, which is not included. 



4 Results and Discussion 

The x° BS distributions shown in Fig. ^ in all three Q 2 ranges, cannot be described 
by HERWIG without including a significant LO resolved photon component, which is 
dominant for x® BS < 0.75. Hence the dijet cross sections in this region are sensitive to 
the photon structure. 

The measured dijet cross sections for the low and the high E^ 1 cuts described in Section 3 
are shown in Figs. ||] and [3|, respectively. The shapes of the dijet cross sections change 
markedly with increasing Q 2 , the cross section in the low-x^ BS region decreasing faster 
than the cross section in the high-x° BS region. This effect is more pronounced for the low 
E 3 ^ 1 cuts. 

The dijet cross sections are compared to the predictions of the HERWIG MC at hadron 



level using different photon PDFs. Those of GRV LO [|J and WHIT2 @] are valid 
for real photons only, have differing gluon distributions and have no suppression of the 
resolved photon component as Q 2 increases. In the SaS ID || model the resolved pho- 
ton consists of two separate contributions, the non-perturbative hadronic 'Vector Meson 
Dominance' component and the anomalous pQCD component, each with different Q 2 
dependence. Specifically, the 'Vector Meson Dominance' component of the resolved pho- 
ton is predicted to decrease approximately as (rn 2 p /(m 2 p + Q 2 )) 2 . The pQCD component 
is predicted to decrease more slowly as ~ log(/z 2 /Q 2 ), where /i 2 is the hard QCD scale 

of the process which, for jet production, is usually taken to be proportional to (E^) . 
The measured cross sections for the LDIS region are also compared to the LEPTO pifl 
Monte Carlo prediction, which does not include a resolved photon component and uses a 
parton-shower model to account for higher-order pQCD effects. The general framework 
is similar to the LO-DIR HERWIG and PYTHIA simulations. In this picture, the dijet 
cross section at low x® BS arises purely from parton-shower contributions to the LO-DIR 
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process. The HERWIG and LEPTO predictions agree in the highest x® BS bin, where 
the direct component dominates. In order to compare the shape of the measured cross 
sections with that of the MC predictions, the latter have been normalized to the data 
cross sections for x° BS > 0.75. 

The low £y e * cross sections are compared to the HERWIG predictions in Fig. |2|. The SaS 
ID prediction without MI agrees qualitatively with the data in the LDIS range; however 
a disagreement is observed at low x° BS in the IQS and PHP ranges, which becomes more 
striking as Q 2 decreases. The GRV prediction without MI and the WHIT2 prediction 
both without and with MI using p™ m = 2.0 GeV are compared to the data in the PHP 
region, where the discrepancy with SaS ID is greatest. The effect of MI is found to be 
very sensitive to the p™ m value and to the choice of PDF |27| , [35|| . The model using 
WHIT2 with MI gives reasonable agreement with the data. The shape of the low E^ e *-cut 
cross sections, shown in Fig. ^| (a), cannot be described by the models that do not include 
MI. The discrepancy seen for the PHP data using SaS ID without MI is not present in 
the LDIS region. This is as expected, in the framework of the MI model, if the resolved 
component decreases with Q 2 . The LEPTO predictions underestimate the dijet cross 
sections at low x® BS in the LDIS region, indicating that the parton-shower contributions 
alone cannot describe the dijet data in this region. 

The high-i?^ 6 * data are shown in Fig. || The predictions of HERWIG without multi- 
parton interactions using the SaS ID photon PDF describe the shape of the measured 
cross section well in the LDIS region but tend to underestimate the PHP and IQS data 
at low x® BS . The measurements are also compared to HERWIG using GRV without MI. 
This model is in good agreement with the data in the PHP and IQS regions but fails to 
describe the data in the LDIS region, as expected since the GRV set describing the real 
photon structure is used. As seen in Fig. |3| (c), LEPTO again underestimates the dijet 
cross sections at low x^ BS . 

The cross-section ratio cr(x° BS < 0.75)/<r(x° BS > 0.75) as a function of Q 2 for both sets of 
Ej? 1 cuts is shown in Fig. f|. The dominant systematic uncertainties of these measurements 
(7-16%) are due to the and x° BS cuts, except for the LDIS samples where the cut on 
the impact point of the scattered positron results in an additional systematic uncertainty 
of about 10%. For the IQS measurements, the latter systematic uncertainty falls to 5%. 
When the data are corrected using PYTHIA, the measured ratios are systematically lower 
for all Q 2 points. This systematic error is therefore not included with the previous ones, 
but it is shown separately. For the PHP data, there is an additional error of 5% due to 
uncertainties in the Monte Carlo normalisation factors for the LO-DIR and LO-RES used 
in the fit (not shown). 

The cross-section ratio falls steeply as a function of Q 2 . This can be interpreted as 
the suppression of the resolved photon component as the photon virtuality increases. 
The decrease is more pronounced for the measurements using the low set of E^ 1 cuts, 
which are more sensitive to the resolved component and a possible underlying event. The 
predictions of HERWIG with two different photon PDFs are also shown. The prediction 
using the GRV set is flat, irrespective of the presence of MI, as expected for a photon 
PDF lacking a Q 2 dependence. The prediction using the SaS ID PDF decreases with Q 2 
and lies below the data in the low Q 2 region. The measured ratios are also compared 
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with the predictions of LEPTO in the LDIS region in which this model is applicable. The 
LEPTO predictions show the contribution to the ratio arising from parton shower effects 
alone and underestimate the measured ratios in both cases. 

In Fig. |(b), the high-E J T et dat a are also compared to the predictions of a NLO pQCD 
calculation at the parton level using the program JetViP The renormalisation and 



factorisation scales were set to Q 2 + (E^ 1 ) . The calculation includes contributions from 
a resolved photon component, which are computed using two different sets of photon 
PDFs: the SaS ID PDF and the GS96 HO |§ PDF modified to include a Q 2 suppression 
according to Drees and Godbole || (GS96 DG). The JetViP predictions are sensitive to 
the choice of the photon PDFs but lie well below the data. The magnitude of the hadron- 
to-parton level corrections has been investigated as a possible source of this discrepancy. 
The data corrections to parton level were estimated using the MC samples and were 
found to decrease the measured cross section ratios by approximately 20-30%, which is 
insufficient to explain the discrepancy. 



5 Conclusions 

Dijet cross sections, da/dx® BS , have been measured using the longitudinally-invariant kx 
jet algorithm as a function of Q 2 , for Q 2 < 1 GeV 2 , 0.1< Q 2 < 0.55 GeV 2 and 1.5< 
Q 2 < 4.5 GeV 2 . The x® BS dependence of the measured dijet cross sections changes with 
increasing Q 2 . The low-x^ BS cross section decreases more rapidly than the high-x° BS 
cross section as Q 2 increases. This effect is more pronounced for the lower of the two sets 
of Ej? cuts. 

The shape of the dijet cross sections, da/dx® BS , is compared to the predictions of HER- 
WIG M C for a variety of photon PDFs. None of these models is able to explain the data 
for both high- and low-i^ 6 * cuts in all Q 2 ranges. 

The ratio a(x° BS < 0.75)/<t(x° bs > 0.75) for dijet cross sections decreases as Q 2 in- 
creases but remains above the level expected from parton-shower effects alone. This 
may be interpreted in terms of a resolved photon component which is suppressed as the 
photon virtuality increases but which remains present up to Q 2 = 4.5 GeV 2 when the 
photon is probed at the scale /i 2 ~ 30 GeV 2 of these measurements. Within the models 
available, events at x® BS < 0.75 can originate from non-perturbative photon structure or 
perturbatively-calculable higher-order processes, and are influenced by underlying-event 
effects especially at low Q 2 and low E 3 ^ . The features and trends seen in the data are 
in accord with general expectations. However, none of the LO models, or the NLO cal- 
culation examined here, gives a good description of the data across the full kinematic 
region. 
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Figure 1: Uncorrected x° BS distributions for dijet events selected with the k T algorithm, 
for E 3 T et > 5.5 GeV and -1.125 < if et < 2.2, in the ranges: (a) Q 2 « GeV 2 , (b) 0.1< 
Q 2 < 0.55 GeV 2 , and (c) 1.5< Q 2 < 4.5 GeV 2 . The points are the measurements and 
the solid histograms are the predictions of the HERWIG Monte Carlo. The simulation 
of multiple parton interactions was used only for the Q 2 ~ GeV 2 sample. In (a) the 
reweighted predictions of HERWIG (dashed histogram) used for the data correction are 
also shown. The shaded histograms represent the LO-DIR contributions. 
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Figure 2: Dijet cross sections, da/dxx. , for jets of hadrons selected with the kx algorithm 
in the Q 2 ranges: (a) Q 2 « GeV 2 , (b) O.K Q 2 < 0.55 GeV 2 , (c) 1.5< Q 2 < 4.5 GeV 2 
for the low E T et set of cuts. The points represent the measured cross sections. The 
inner error bars represent the statistical errors, and the outer are the statistical and 
uncorrelated systematic errors added in quadrature. The shaded band represents the 
systematic uncertainty due to the modelling of the jet fragmentation, estimated using 
PYTHIA. The shaded horizontal band represents the uncertainty due to the CAL energy 
scale. The full histogram represents the HERWIG predictions without MI using the SaS 
ID photon PDFs, and the dotted histogram represents those with the GRV LO real 
photon PDFs. The predictions of HERWIG without MI using the WHIT2 set (dashed 
histogram) and with MI for p™ in = 2.0 GeV (dot-dashed histogram) are shown in (a). 
The predictions of LEPTO are shown in (c) as the dot-dashed histogram. The LEPTO 
and HERWIG predictions have been normalised to the x® BS > 0.75 dijet cross section. 
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Figure 3: Dijet cross sections, da/dx® , for jets of hadrons selected with the kx algorithm 
in the Q 2 ranges : (a) Q 2 « GeV 2 , (b) O.K Q 2 < 0.55 GeV 2 , (c) 1.5< Q 2 < 4.5 GeV 2 
for the high E T et set of cuts. The points represent the measured cross sections. The 
inner error bars represent the statistical errors and the outer are the statistical and the 
uncorrelated systematic errors added in quadrature. The shaded band represents the 
systematic uncertainty due to the modelling of the jet fragmentation, estimated using 
PYTHIA. The shaded horizontal band represents the uncertainty due to the CAL energy 
scale. The full histogram represents the HERWIG prediction without MI using the SaS 
ID photon PDFs, and the dotted histogram the GRV LO photon PDFs. The predictions 
of LEPTO are shown in (c) as the dot-dashed histogram. The LEPTO and HERWIG 
predictions have been normalised to the x® BS > 0.75 dijet cross section. 
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Figure 4: The ratio of dijet cross sections, cr(x° BS < 0.75)/cr(x° BS > 0.75), as a function 
of photon virtuality, Q 2 , for dijet events selected using the kx algorithm. The points 
represent the ZEUS data. The cross section ratios are shown for both the low (a) and 
the high (b) E 3 ^ 1 cuts. The inner error bars represent the statistical errors and the 
outer are the statistical and systematic errors added in quadrature. The shaded band 
represents the systematic error due to the uncertainty in modelling the jet fragmentation, 
estimated using PYTHIA. The shaded horizontal band represents the uncertainty due 
to the CAL energy scale. Also shown are the predictions of HERWIG without MI for 
two different choices of photon PDFs: GRV for real photons (dotted histogram) and SaS 
ID (full histogram). The LEPTO predictions are shown for Q 2 > 1.5 GeV 2 (dot-dashed 
histogram). The predictions of NLO pQCD calculated using the JetViP program with 
SaS ID Photon PDFs are shown as the full histogram and those using GS96 DG are 
shown as dashed the histogram. 



15 



